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Abstract. We study convexity and monotonicity properties for prices of bonds and 
bond options when the short rate is modeled by a diffusion process. We provide 
conditions under which convexity of the price in the short rate is guaranteed. Under 
these conditions the price is decreasing in the drift and increasing in the volatihty of 
the short rate. We also study convexity properties of the logarithm of the price. 



Already in the seminal paper convexity of the option price in the underlying 
asset is discussed. In the last decade this issue has attracted renewed interest in the 



literature, compare [T], [1], [3], 0, [H], [13], [H], [E], [H!, [IE], [H] and [20]. Given 



a convex pay-off function one asks for what models this convexity is preserved in the 
sense that the price also is a convex function of the underlying asset at any fixed time 
prior to maturity. This question is studied for diffusion models, for models with jumps 
and for various option types including options written on several underlying assets. 
The interest in convexity has at least three reasons. Firstly, convexity is a fundamental 
qualitative property of option prices (an even more fundamental qualitative property 
would be monotonicity in the underlying asset provided the pay-off is monotone, but 
such properties can usually be derived immediately). Secondly, if the price is convex 
then it is also typically increasing in the volatility, and in the case of jump-diffusion 
models, also in the jump parameters. Thirdly, if a delta-hedger uses a model that over- 
estimates the true volatility, he or she will obtain a superhedge for the claim provided 
the price is convex. 

Our aim with the present paper is to continue this study to bonds and bond options, 
for which we regard the short rate as the underlying process. Thus we study preserva- 
tion of convexity for the term structure equation instead of variants of the Black-Scholes 
equation as is the case in the references above. Surprisingly little has been done in this 
direction, with [2] as a notable exception. A motivation for this is perhaps that the 
third reason for studying convexity mentioned above is not directly applicable since the 
short rate is not a traded asset. However, the first two stated motives remain valid also 
in interest rate theory. 

We assume that the short rate is modeled under some given risk neutral probability 
measure as a stochastic process X = (Xt)t>o with dynamics 



where B is a standard Brownian motion and /3 and a are given functions of time and 
the current short rate. We first investigate the convexity properties of the T-bond price 



1. Introduction 



dXt = /3{Xt,t)dt + a{Xt,t)dB, 
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where the indices indicate that Xt = x. Using the Feynman-Kac theorem it follows 
that the bond price u satisfies the term structure equation 



This enables us to use a mixture of stochastic techniques and methods from the theory 
of parabolic partial differential equations. 

One should note that our pay-off function in the case of bonds is identically equal to 
one, so it is both convex and concave. However, convexity in x is the natural property 
to consider: since the bond price declines with x, convexity means that the absolute 
value of this decline decreases with x. We certainly expect bond prices to suffer a 
smaller decline if short rates move from 5% to 6% than if the short rates move from 
1% to 2%! In Section [S] we find conditions on the model for X that guarantee that 
convexity is preserved. More precisely, if (3xx ^ 2 (in the sense of distributions if /? 
is not twice continuously differentiable) , then the model is convexity preserving. To 
our knowledge, this condition is indeed fulfilled for all models of the short rate that 
are used in practice, compare Table [1] below. Moreover, the condition is sharp in the 
sense that it is also a necessary condition for preservation of convexity provided the 
coefficients of the model are regular enough, see Theorem l5.2[ Using the pathwise non- 
crossing property of diffusions, it is easily seen that the bond price u is decreasing in 
the drift /3. In Section [6] we show that the bond price in a convexity preserving model 
is also increasing in the volatility a. Thus the general relationship between convexity 
and monotonicity in the volatility known from option pricing theory extends to interest 
rate theory. 

Our study of convexity and monotonicity properties is formulated for prices of options 
written with the short rate as the underlying asset. Thus we study the function 



where g is a given convex pay-off function (note that the case g = 1 corresponds to a 
bond). The main reason for extending the study from bonds to options on the short 
rate is to be able to study bond option prices, i.e. prices of options written on a bond 
price as the underlying asset. In fact, our general convexity and monotonicity results 
allow us in Section [7| to deduce properties of certain bond options. For example, the 
price of a bond call option is convex in x and therefore also increasing in a for convexity 
preserving models. 

It is also natural to consider convexity properties of the logarithm of the bond price. 
This is connected to the notion of duration, i.e. the negative of the derivative of the 
logarithm. The analogous concept for stock options is elasticity, compare [S], [TU], [TC] 
and [22] . We say that the price is log-convex if the logarithm of the price is convex in x 
and analogously for log-concavity. Again, since the pay-off is constant for the bond it is 
both log-convex and log-concave. Unlike the case of convexity, however, both of these 
cases deserve consideration. According to the discussion above, convexity of a bond 
price means that the absolute value of the decline is decreasing in x. In contrast to 
this, log-convexity means that the relative decline diminishes when x grows (a declining 
duration), and log-concavity means that the relative decline of the price increases (an 
increasing duration). In Section [8] we show that if the drift f3 is spatially concave and the 
diffusion coefficient o"^ is spatially convex, then log-convexity is preserved. Similarly, in 
Section[9]we show that log-concavity is preserved provided f3 is convex and o"^ is concave. 
If we insist that the model should preserve both log-convexity and log-concavity, we 
arrive at models where the logarithm of the bond price is both convex and concave. 




for t < T 
for t = T. 
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i.e. linear. These are, of course, the models that admit an affine term structure (in 
our context these bond prices would be referred to as being log-affine). Apart from 
admitting explicit bond prices, affine models play an important role in interest rate 
theory, compare for example |8i]. Thus we recover the well-known sufficient condition 
that P and o"^ are affine for the existence of an affine term structure. 

In the next section we present the assumptions on the model parameters under which 
our results are presented. If additional regularity of the coefficients is assumed, bounds 
on the spatial derivatives of bond and option prices can be obtained, see Section [3l In 
Section H] a continuity result is provided, thus allowing us to assume that the coeffi- 
cients are regular, and the general results follow by approximating the coefficients. As 
mentioned above, our main results are presented in Sections [MHl and a summary of 
these results is presented in Section [lOl 

One important feature of our approach is that it works just as well for models where 
X never reaches zero, for models where the rate can reach zero, as for models that 
allow negative interest rates. We thus avoid a case by case study. However, we will 
throughout the paper point out to which of the commonly studied models presented in 
Table [Hour results are applicable. For a more detailed discussion of these models, see 
for example 



Model 


Dynamics 


X>0 


AB 


AO 


V 


dX = k{d - X)dt + adB 


No 


Yes 


Yes 


CIR 


dX = k{e - X)dt + aVX dB 


Yes 


Yes 


Yes 


D 


dX = hXdt + aX dB 


Yes 


Yes 


No 


EV 


dX = X{r] - a\nX) dt + aX dB 


Yes 


No 


No 


HW 


dX = k{et - X) dt + a dB 


No 


Yes 


Yes 


BK 


dX = X{r]t-alnX)dt + aX dB 


Yes 


No 


No 


MM 


dX = X{r)t - (A - In X) dt + aX dB 


Yes 


No 


No 


Table 1. Some examples of short rate models. ^- 


?he table is copied 



from [llj. AB stands for analytic bond prices, and AO stands for ana- 
lytic bond option prices. When referring to these models below we will 
assume that k, 6, a, rj and 7 are positive and that A > 7. 



2. Assumptions 

As explained in the introduction, we assume that the short rate X is modeled as 

(3) dXt = (3{Xt,t)dt + a{Xt,t)dBt 

for some Brownian motion B. Note that some of the short rate models in Tabled] are 
specified so that the interest rate cannot fall below zero, whereas some models allow 
for negative interest rates with positive probability. To unify the analysis, we view the 
short rate process X in ([3]) as specified on the whole real line M with a and /3 suitably 
extended (for example to be 0) for negative values in case X is specified initially only 
on the positive real axis. 

Throughout the article we make the following regularity and growth assumptions. 

Assumption 2.1. The functions /?, cr : M x [0, cxd) M are continuous functions, (3 
is locally Lipschitz continuous in the x-variable, and a is locally Hdlder(l/2) in the 
x-variable. Moreover, there exists a constant D such that 

(4) W{x,t)\ < D{l + x^) 
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and 

(5) |/3(x,t)| <Z)(l + |x|) 

for all {x,t) G R X [0,00). 

The conditions on a and (3 guarantee a non-exploding unique strong solution to 
([3]) for any initial point {x,t) G M x [0, 00). The condition that a is bounded for 
negative x implies that, for any T, the price u in ([T]) of a T-bond is finite at all points 
{x,t) £ M X [0,T], see Corollary 13.31 below. Without this condition, the bond price 
u may be infinite. Indeed, models in which a grows faster than ^/\x\ for negative x 
have typically infinite bond prices, compare Theorem 4.1 in [26]. One should note 
that all models in Table [T] satisfy the conditions of Assumption 12.11 in the following 
sense: those models in which the short rate may fall below zero (V and HW) clearly 
satisfy Assumption I2.lt the remaining models give rise to non-negative rates, so the 
parameters can be chosen arbitrarily for negative values so that Assumption 12.11 holds 
without affecting the bond value for positive rates. 

3. Auxiliary estimates 

When dealing with bonds and options on the short rate it is natural to study the 
corresponding term structure equation, i.e. the terminal value problem 

Ut + aU^^ + pUa: -xU = for t < T 
U = g for t = T, 

where a := (t^/2. Due to technical reasons we instead study the function 
(7) V{x,t):=Vf{x,t):=E,^t 



(6) 



expj-^ fiXs)ds^g{XT] 



(8) 



for some appropriate function / and then let / approach x. Note that the corresponding 
terminal value problem is 

Vt + aV^, + pV,-fV = for t < T 
V = g for t = T. 

Our choice of / is indicated in the following hypothesis. 

Hypothesis 3.1. / is smooth, concave and there exists a constant K' > such that 

X if X < K' 

constant if x > K' + \. 

The main goal of this section is to provide some estimates on option prices and the 
derivatives of V , compare Corollary 13.31 and Proposition 13.51 We first claim that V 
grows at most exponentially as x — > —00. To be more precise, assume that the pay-off 
function g satisfies 

(9) < g{x) < Mmaxje"-^"^, 1} 

for some non-negative constants M and K (if g = 1 then ([9]) holds with M = 1, K = 0). 
Moreover, let D be the constant from Hypothesis 13.41 From ([5]) we have 

(10) P{x, t) > —D{1 — x) for all negative x. 
Define 

(11) M/(x,i) :=e^(^')^(*V(x,t), 
where 



eD{T-t) _ 1 



h{t) = - _^ + Ke^(^-*). 
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The choice of the function h is motivated by the following observation. If (t(x, t) = 
and (5{x,t) = Dx for negative x, and if Xq = x < 0, then X satisfies Xt = xe^^ . 
Consequently, 



V{x,t) 



E. 



x,t 



< exp < —X- 



exp |- fiXs) ds^ giXr) = {' ^^''^"'^ ds^ gixe""^^-'^] 

D{T-t) _ I ' 



D 



M 



exp |-Kxe^(^"*)| = Me" 



■xh{t) 



where we have used the bound ([9j). This indicates that the function W is bounded (at 
least for negative values of x). As is shown in Lemma [32] below, this intuition is indeed 
correct. 

Lemma 3.2. Assume Hypothesis \ 3. 1\ and the bound ([9]). Then the function W defined 
in (jlip is bounded on M x [0,T], i.e. there exists a constant C such that <W < C . 

Proof. By the bound ([9]) and the Cauchy-Schwartz inequality we have 
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W^ixA) 



Jix)h(t) 



M2 \ M 

(12) < (^E,^t 
< Ex,t 



-Ex,t 



expj-^ f{Xs)ds^g{XT) 



exp 



Ms-t) j(^) _ f^x,) ds \ max{e-^^-, l}e^/(-')^"*"'*' 



exp|2y|'^e^(^-*)/(x)-/(X,)ds 
max{e-2^^-,l}e2^^(^)^"*""'^" 



By Jensen's inequality (applied to the exponential function and the time integral) we 
have 

rT 



E. 



x,t 



exp|2^ e^(^-*)/(x)-/(X,)ds 

exp {2(e^(^-*)/(x) - /(X,))(T - t)] ds 



< Ex,t 
1 



T 



E 



exp \ —2e 



D(s-t) 



{T-t)Y,] 



ds 



T-tJt 

where Yg = e~^^^^^^ f{Xs) — f{x) with 1^ = 0. By Ito's lemma we have 

dYs = Psds + ds dBs for s >t, 

where 



Ps = e-^(^'*) ( f'iXsmXs,s) - Df{Xs) + lnXs,s)a\Xs,s) 



and 



o-.o = e 



f'iXs)aiXs,s) 



It follows from the assumption (jlOp and Hypothesis [3TT] that /3 is bounded from below, 
i.e. there exists a constant C such that (ig^ —C for all s almost surely. Similarly, the 
bound dH) and Hvpothesis 13.11 vield that \as\ < C. Since y ^ exp{-2e^('^-*)(T-t)y} 
is decreasing and convex, we have 



(13) 



E 



exp{-2e^(^-*)(r-t)n} <E exp{-2e^^^-'^(r - t)y,} 



oD{s-t), 
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for every s £ [t,T], where 

= -C{s -t) + CBs^t 

is a Brownian motion with drift starting from at time t. Indeed, monotonicity in 
the drift is immediate, and monotonicity in the volatihty for processes with constant 
drift and convex pay-off functions is weh-known, compare Theorem 6.2 in [T^ and 
Theorem 7 in [20]. Consequently, 



E. 



x,t 



exp|2^ /(x)e^(^-*)-/(X,)ds 



< 



-l^"" E [exp {-2e^(^-*)(r - t)Y,] 



ds < C 



for some constant C (in fact, Yg is N{—C{s — t),C'^[s — t))-distributed, so exphcit 
bounds of the above integral are readily derived). Consequently, the first factor in (fT2]l 
is bounded. For the second factor in (fT2]l. note that 



(14) 



'x,t 



maxje 



-2KXt 



l}e 



2Kf{x)e°^T-t) 



< E. 



x,t 



,2K(f(x)eO(T-'^^f{XT)) 



+e 



The first term in ()14p is of the form shown above to be bounded, and the second one is 
bounded since / is bounded from above. This shows that W is bounded, which finishes 
the proof. □ 

Corollary 3.3. Assume that g satisfies ([9]). Then the option price U defined in ([2]) 
above is finite. In fact, there exists a constant M > such that 

U{x,t) < Mmax{l,e~^'^''} 

for all {x,t) e M X [0,r]. 

Proof. Since f{x) < x, we have that U{x,t) < V{x,t) for all x and t. Consequently, it 
follows from Lemma 13.21 that for every choice of a function / satisfying Hypothesis 13.11 
there exists a constant C such that U{x,t) < Ce~^^^^'^^^'^'^ *)-i)/D-i^e^(^ for all 

(x,i) G M X [0,r]. Thus 

U{x,t) < Mmax{l,e-*^''} 
for some large constant M. □ 

We next provide conditions under which the A;th spatial derivative of the function 
W decays like jxl"*^, k = 1,2,3. We need the following assumptions. 

Hypothesis 3.4. The coefficients a and P are smooth and a > at all points. More- 
over, 

(15) \P{x,t) — Dx\ does not depend on x for {x,t) G M x [0,T] with x large negative, 
and the derivatives satisfy the growth conditions 

|a^'/3(x,t)| <C(l + |x|)i-^ 
for {x,t) G [0,cx)) X [0,r], 

\d^a{x,t)\ < C7(l + Ixl)^-'^ 
for {x,t) G [0,oo) X [0,r], and 

\d^aix,t)\ < C7(l + Ixl)-'^ 
for {x,t) G (-00,0] X [0,r]. 



(A: = 0,1, 2, 3) 
(fe = 0,l,2,3) 
{k = 0,1,2,3) 



CONVEXITY THEORY FOR THE TERM STRUCTURE EQUATION 



Proposition 3.5. Assume Hypotheses \3.1\ and \3.4\ that the pay-off function g is 
smooth, satisfies Q and that e^^^^^ g{x) is constant for large \x\. Then there exists a 
constant C such that the function W defined in satisfies 



(16) 

for k = 0,1,2,3. 



\d^w\ < ,. 



Proof. The case A; = follows from Lemma 13.21 Thus we know that 1^ is a bounded 
solution to 



(17) 
where 

and 



Wix,T)=gix) 

5(x) = e^(^)^5(x), 
P = I3- 2fi,ah 



l={Df- f^l3)h + f'.ah' - f^^ah. 
Let W{x,t) = W{x,t) -g{x). Then 



(18) 
where 



W{x,T) = 0, 



h = ag^x + Pgx + 75- 

Since g is constant for large values of W satisfies if and only if W does. Using 
Hypotheses 13.11 and 13.41 we note that 

\d^a\ < C{1 + \x\f-'' \d^P\ < C(l + |x|)i-*^ 
\d^7\ < C(l + |x|)"*^ \d^h\ < C(l + \x\)-'' 

for some constant C and for all k = 0,1,2,3. By the Feynman-Kac representation 
theorem, 



W{x, i)=Ey expJ^j 7(y/='*, r) dr j /i(y/'*) ds 
where Y^'^ is the diffusion 

'* = /3(y/'*, s) ds + a(Yf'\ s) dBs 
with Y^'^ = X. Now, let x < y. Then 



\W{x,t)-W{y,t)\ < E 

cT 



exp 



7(y;^'*, r) dr } h{Yf'') - h{Yy'') 



+E 



exp 



7(y,:"'*, r)dr \ - exp <^ / 7(^/'*, r) dr 



ds 

h{Yy'^) ds 



= h + h. 

Since 7 is bounded and h is Lipschitz continuous, we find that 

h<c r E [|y/'* - y/'*|] ds = c r E [Yf - yf '*] ds, 

Jt Jt 



8 



ERIK EKSTROM AND JOHAN TYSK 



where the equahty holds since x < y implies Yf'^ < Ys'^ for s > t, compare Theo- 
rem IX. 3. 7 in [23] . Since the drift /? is Lipschitz continuous, we have 

E [Yy'' - y/'*] = y-x+l' E r) - P{Y,-'\ r)] dr < y-x+C j^^ E ^^^^ - y,-'*] 

for any s > t. It follows from Gronwall's lemma that 

Ii<C\y-x\{T-t). 

Similarly, since h is bounded and 7 is bounded and Lipschitz it can also be shown that 

h < C\y - x\{T - t) 
for some constant C. It follows that Wx is bounded on M x [0, T] and that 

lim Wx(x,t) = 0. 

{x,t)^{xo,T) 



(19) 



dr 



Therefore, by differentiating the equation (jlSp we know that the first spatial derivative 
p := Wx is a bounded classical solution to 

Pt + apxx + (ax + 0)Px + (7 + Px)p + IxW -hx = 
p{x, T) = 

It is straightforward to check that 

p{x,t) := e^^^-'h{x), 

where I is a smooth positive function with l{x) = 1/(1 + |x|) for large is a super- 
solution to ()19p for some large constant C, and —p is a subsolution. By the maximum 
principle we have \Wx\ = |p| < P < e'^"^/(l + |x|). This proves (fT6]) for k = 1. 

Next, the solution to equation (I19p also has a stochastic representation. Using this 
representation, the same analysis as above can be applied to show that px is bounded, 
and using the maximum principle again it is straightforward to check that Px = Wxx 
decays like Finally, the same reasoning can be applied to prove (I16p for k = 3. □ 



4. Continuity of bond prices in the model parameters 

As mentioned in the introduction, it is central to our approach to be able to approx- 
imate a bond price (or an option price) by a sequence of bond prices (option prices) 
in models with regular coefficients. Thus we need to know that bond prices are con- 
tinuous in the model parameters. To formulate such a result, let (T„ and /?„ satisfy 
the Assumption 12.11 uniformly in n (i.e. the bounds (jH) and ([5]) hold with the same 
constant D). Also assume that — ^ o and Pn ^ P uniformly on compacts as n — > 00. 
Let X" be the solution of the stochastic differential equation 

dX^ = Pn{X^, s) ds + an{X2, s) dBs X]^ = x. 
Under these conditions it is known that 



(20) Ex^t sup 

t<s<T 

as n ^ 00, see Theorem 2.5 in |3J. Now let 

(21) Un{x,t) = Ex,t' 







expj-^ X^ds^g{XJ}^] 



Here g is assumed to have the following property: for each positive constant k there 
exists a constant such that 

(22) \g{x)^k-g{y)^k\<Ck\x-y\ 
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for all x,y € M. Note that this condition is satisfied for example if g is convex and 
decreasing. 

Proposition 4.1. Assume that an, Pn <ind X" are as described above. Also assume 
that g satisfies ^ and (|22|). Then 



U{x,t) = lim Un{x,t) 



for all {x,t) E M X [0,T]. 



Proof. Fix a point G M x [0,T]. For positive constants k, define U^{x,t) and 

U^{x,t) by 

rT 

Tk, 



and 



Then 



where 



and 



U''{x,t) = E^^t expj^ kA-X,ds^{g{XT)Ak) 

U^{x, t) = E^^t exp {^^ A -X; dsj (5(X^) A k) 
U^{x,t)-U\x,t) 
exp!. [ kA -X^ ds\ \giX^T) A ^ - giXr) A k\ 



<h + h, 



1 1 — Ex 



h — Ex,t 



expi / kA-X^dt 



k A -X. ds 



(giXr) A k) 



Note that 

h < e^^^-'^'E^^t [|5(^t) a A; - giXr) A k\] < e^^'^-'^CuE^^t [|X^ - Xt\] 

since g satisfies (I22p . It thus follows from (I20p that /i ^ as n — > cxd. Similarly, using 
the inequality |e^ — e^| < (e^ + e'^)\x — y|, we find 



h < kE^^t 

< 2ke^^^-^^ E^^t 

< 2ke^^^-'^E,^t 
as n — > oo by ([20|) . Consequently. 



expi^J^ /c A -X" dsj - exp 1^ fcA-Xgdsj 



T 



|A:A-Xr-A;A-XJ(is 



(T-t) sup |X,"-X, 

t<s<T 



(23) lim = U^{x,t) 

n^oo 

for each fixed k. It then follows from the monotone convergence theorem that 

(24) lim lim [/^(x,t) = lim E^^t 



k—*oc n^oo 



Since 



k^oo 
Ex,t 



Un{x,t) = lim U^{x,t) 

fc— >oo 



expj^ kA-f{Xs)dsY9iXT)/\k) 
exp!^- j^^ fiX,)ds^giXT) =U{x,t). 
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by the monotone convergence theorem, it suffices to show that the order of the hmits 
in ()24p can be interchanged. To see this, first note that 

(25) < Un{x,t)-U^{x,t) 



+Ex,t 

We first claim that 
(26) 

To see this, note that 
< Exf 



exp |- X2 ds^ giX^) - exp A: A -X^ dsj giX^) 
exp |- X^ dsj ig{X^) - k A giX^)) 



hm sup li = 



h + h- 



for i = 3, 4. 









exp < 


















exp < 













< 



E. 



x,t 



exp 



X^dsUg{X^)f 



Ex,t 



j\-k-X^)^ds^ ^ 



1/2 



where we have used the inequahty 1— e^' < —x and the Cauchy-Schwartz inequahty. The 
first factor is bounded uniformly in n for each fixed x (this can be proved analogously 
to Lemma l3.2p . As for the second factor, an application of Jensen's inequality gives 

ds. 



Ex,t 



{-k - X^)+ ds 



<{T-t) Ex,t[{{-k-X^)+)' 
Proceeding as in the proof of Lemma 13.21 above, it is straightforward to check that 



^ 2DT p 

S e ±Jx,t 



for a Brownian motion Y with negative drift. Moreover, since {{—ke — Ys)'^)'^ is a 
submartingale we have 



Ex,t 



i-ke 



-DT 



Y' 



<Ex-, 



i-ke 



-DT 



YtY 







as A; ^ oo by dominated convergence. Consequently, (|26p holds for i = 3. Similar 
methods can be employed to establish (j26p also for i = 4. 

Let e > be given. From ([25]) and (f26l) it follows that there exists a ko such that 

< Un{x,t) - U^{x,t) < e 

for all k > kQ and all n. Thus, for such a k we have 

lim U^{x,t) < liminf C/„(x,t) < limsup C/„(x, t) < lim Un{x,t) +e 

n-^oo n^oo n— >oo n— +00 

where we recall from ()23p that the outer limits exist. Letting A; — > cxd and remembering 
(j24p above we find that 

U (x, t) < liminf Un{x, t) < lim sup Un{x, t) < U (x, t) + e. 
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Since e is arbitrary, 

lim Un{x,t) = U{x,t) 

n— >oo 

as required. □ 

5. Convexity of bond prices 

In this section we show that interest rate models have convex bond prices provided 
the drift /? is not "too convex". More precisely, we need to require that Pxx < 2 (in 
the sense of distributions). To see why this condition comes into play, consider the 
corresponding term structure equation, i.e. the terminal value problem 

Ut + aU^x +pUx-xU = 
Uix,T)=g{x) 

for some convex and decreasing pay-off function g. When using the PDE- approach we 
will in the sequel simplify the presentation by performing a standard change t = T — t 
of the direction of time. By a slight abuse of notation, we use the same symbols /?, a, 
U, V etc. to denote the new functions depending on time r to maturity rather than on 
the actual time t. With this new convention, the term structure equation becomes an 
initial value problem 

' Ur = aUxx + PUx - xU 
Uix,0)=gix). 

Assume for the moment that all coefficients are regular enough. Also assume that there 
is a first point (xo,to) at which convexity is almost lost, i.e. that U{x,t) is convex for 
< r < To and Uxx{xo,tq) = 0. Then 

drUxx = dlUr = dliaUxx + PUx - xU) 

= aUxxxx + (2aa; + P)Uxxx + (2/3x - x)Uxx + (Pxx - 2)C/x- 

Since x i-^ Uxx{x,to) has a minimum at x = xq, we llQjVG ^-^xxxx ^-^xxx — ^-^xx — ^ Sit 
(xq, To). Thus we find that 

drUxx > {Pxx - 2)Ux > 

at this point provided Pxx ^ 2, i.e. the infinitesimal change of U at the critical point 
(xo,To) is convex. This suggests that convexity will be preserved if Pxx < 2. Below we 
make this argument rigorous. 

Theorem 5.1. Assume that Pxx{x,t) < 2 (in the sense of distributions) at all points 
{x,t) £ M X [0, T]. Also assume that the pay- ojf function g is convex, decreasing and 
satisfies ([9]). Then the corresponding option price U{x,t) is convex in x at all times 
t£[0,T]. 

Remark Note that all models in Table [J satisfy the condition Pxx < 2. Consequently, 
all those models give rise to convex bond prices. 

Proof. We first assume that Hypothesis 13.41 holds, that Pxx = for x > C for some 
constant C, and that the pay-off function g is smooth with e^-^^^^ g{x) being constant 
for large |x|, where K is the constant from ([9]). 

Instead of studying the option price U directly, we first study the function V = 
defined in ([7]) for some function / satisfying Hypothesis 13.11 such that 

(27) l3xx{x,T)-2fx{x)<{) 

for all (x,r) G M x [0, T]. Under these assumptions we claim that the function V is 
spatially convex. 
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To see this, note that it follows from Proposition 13.51 that there exists a constant C 
such that 

(28) \V,,{x,t)\<Cp{x,t) 
at all points, where 

(here Kq is a large constant so that Kq > (e^'^ — l)/D + Ke^^). For e > 0, consider 
the function 

V%x,t) := V{x,T) + ee^^^x'^ + x^p{x,t)). 

Here the even number > 2 is chosen large so that + x^p{x,t) has a strictly 
positive second spatial derivative at all points (x, r) G M x [0, T] (the constant M will 
be chosen below). Since V satisfies drV = C^V where 

£/ = adl + /3d^ - f, 

we find that 

(29) dlidrV - C^V) = edl{dr - Cf) {e^^^{x^ + x^p)) = ee^'^^{Mh - h), 
where 

Ii = dlix'^ + x^p) > 

and 

h = dl (adl{x^ + x^p) + m^{x? + x^p) - fix? + x^p) + (£> + \)K^e^^^^^^fx^' 

For large positive x we have that / is constant, so p is bounded. Thus the term I\ grows 
like whereas the term I2 grows at most like x^~'^ for large positive x. Similarly, 

for large negative x we have f[x) = x, so 

h ~ Ky^^'^^^^x^'p, 

whereas 

h ~ 0x^0=0=0= + (-/3i^oe^''+'^" - X + (D + l)i^oe(^+')"x)x^fc. 
Using ([5]) we find that the highest order terms of I2 behaves at most like 
+ (-DKoe(^+i)^ - 1 + (D + l)/^oe(^+')^)x^+W, 

and since —DK^e^^^^^'^ — 1 + (D + l)Koe^^^^'*^ > we find that there exists a positive 
constant C such that 

h < Cx^p 

for large negative x. Consequently, Ii grows at least as fast as I2 for large values of 
so M can be chosen large so that M/i dominates I2 everywhere. Actually, we will 
assume that M is chosen so large that 

(30) M/i >h-h 
at all points, where 

/3 := (/3xx- 2/^)5, (x^ + x^p) 
(compare (|32p below). This can be done since f3xx = and thus Is = for large x, and 
/a > for large negative x, compare ([23) ■ 
Now, let 

r:={(x,r) :y/,(x,r) <0}. 

We claim that the set F is empty. To see this, note that the second spatial derivative of 
x^ + x^p grows at least like Cx^p as |x| — > 00. Consequently, it follows from (j28p that 
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r C [—R,R] X [0,r] for some R > 0. Thus F is bounded and T is compact. Suppose 
that r 7^ 0, and define 

To := inf{r > : (x, r) G F for some x € M}. 

Due to compactness, the infimum is attained at (xo,to) for some xq, and by continuity 
of Vxx we have ^^^.(xcro) = 0, so tq > 0. Since V^^.{xo,t) > for r < tq, we must 
have 

(31) dldrV%Xo,To) = drdlV'{xo,To) < 0. 

Moreover, straightforward calculations give 

dliL^Vn = dl{aV^, + l3V^-fV') 

Hi^xx - - fxxV. 

Since V^^{xQ,tQ) = and is convex, the function x i— > V^^{x,tq) has a minimum 
point at Xq. Consequently, V^^^{xo,tq) = and V^^^^{xo,tq) > 0. Therefore, at the 
point (xqjTq) we have 

(32) dliCfV') = aV^,,, + {(3xx-V.)V^-UxV' 

> ee'^^flxx - 2fx)dx{x^ + x^p) = ee^'^h 

where we have used Vx < 0, (f271l and that / is concave. Combining (f29l) . (f30]) . (fSTI) 
and ([32]) yields 

a2(a,y^ - £-^y") = ee^^Mh - h) > -ee^^h > dl{drV' - C^V) 

at (xo,To). This contradiction shows that the set E is empty, i.e. is convex. Letting 
e I it follows that also V is spatially convex at all times r € [0,T]. 

To deduce the convexity of C/, consider an increasing sequence of functions 

fi satisfying Hvpothesis 13.11 such that /j(x) — > x as i ^ oo for all x. Then 

/ h{Xs)ds^ I Xsds 
Jt Jt 



almost surely as i — > oo. Therefore 

Vf'{x,t) = E, 



■^x,t 



x,t 



expj-^ n{Xs)ds'^g{XT 
expj-^ Xsds^giXr) 



U{x,t) 



as i — > oo by monotone convergence. Since each function V-^^ is spatially convex it 
follows that the option price U{x,t) is spatially convex. 

To remove the Hypothesis 13.41 and the assumption that Pxx = for large x, assume 
only Assumption 12.11 and that Pxx ^ 2 (in the sense of distributions). Let {an}'^=i 
and {Pn}^=i be sequences of smooth continuous coefficients satisfying Hypothesis 13.4 
such that {Pn)xx = for large x and {f3n)xx < 2. Moreover, assume that cr„ and /3„ 
satisfy the growth conditions ([!]) and ([5]) uniformly in n and that an ^ cr and f3n ^ P 
uniformly on compacts as n ^ oo. It follows from Proposition 14.11 that 

U{x,t) = lim Un{x,t), 



where the function Un is defined as in (I2ip . Since the pointwise limit of a sequence of 
convex functions is convex, it follows that U is convex. 

Finally, to remove the assumptions about the smoothness of g and that ^^^^ g{x) 
is constant outside a compact we approximate g from above by a sequence {^nj^i 
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of smooth pay-ofFs behaving hke e~^^^^^ outside compacts, such that gn{x) [ g{x) 
as n ^ oo. By monotone convergence it foUows that U is convex also without the 
smoothness requirements on g. □ 

The heuristic calculations in the beginning of this section indicate that the condition 
Pxx < 2 is not only a sufficient condition, but also a necessary condition for preservation 
of convexity for the term structure equation. Our next result shows that this is indeed 
true provided the coefficients are regular enough. 

Theorem 5.2. Assume that a and [5 are smooth. Also assume that a > and Pxx > 2 
at some point (xq , T) . Then there exists an option with maturity T and with a decreasing 
and convex pay-off g such that the corresponding price U{x, t) is non-convex at some 
time t < T. 

Proof. Let g be a smooth convex pay-off function which is linear and strictly decreasing 
in a neighborhood of xq. Since U{x,t) is a solution of a parabolic differential equation 
with regular coefficients, its derivatives exist and are continuous up to the boundary 
r = 0, see [23] (here we again let t = T — t). Straightforward calculations yield 

drU^x = dlUr = dliaU^x + f3Ux - xU) 

— O/Uxxxx 

+ (2a^ + l3)Uxxx + {a^x + WxW^x + (Pxx - 2)^7^. 
Since g is linear in a neighborhood of xq, we find that 

drUxx{xo,0) = {Pxx - 2)fifa' > 0. 

Since Uxx{xq, 0) = 0, this means that U is not convex at some time r > 0. This finishes 
the proof. □ 

6. Parameter monotonicity 

In this section we utilize the well-known connection between convexity and parameter 
monotonicity, see for instance [I^ or [14J . We thus show how preservation of convexity 
implies that bond prices and prices of convex options are monotonically increasing in 
the volatility. To formulate the result, let X and X be two diffusion processes satisfying 

dXt = (3iXt,t)dt + a{Xt,t) dBt 

and 

dXt = P{Xt,t) dt + a{Xt, t) dBt, 

respectively. Let 



and 



U{x 


t) 


— Ex^t 


exp < 


[-fx..,,] 


■ aiXr) 


U{x 


t) 


= Ex^t 


exp < 




■ aiXr) 



be the corresponding option prices. 

Theorem 6.1. Assume that P{x,t) < P{x,t) and \a{x,t)\ > |fT(x,t)| for all {x,t) £ 
M X [0, T], and that either Pxx < 2 at all points or Pxx < 2 at all points (both in the 
distributional sense). Also assume that the pay-off function g is convex, decreasing and 
satisfies Then U{x,t) < U{x,t) for all {x,t) E M x [0,r]. 
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Remark Monotonicity in the drift is immediate since g is decreasing and it is known 
that (5{x,t) < (5{x,t) for all x and t implies that Xt < Xt (if a{x,t) = a{x,t) for all 
X and t). Theorem 16.11 tells us that bond prices are also increasing in the volatility, 
provided the condition f3xx < 2 is fulfilled. In particular, all models in Table [1] are 
monotonically increasing in the volatility and decreasing in the drift. 

Proof. Let / be a given function satisfying Hypothesis 13. H and define the functions 
V{x,t) = Ex,t' 



and 

V{x,t)=Ex„ 



expj-^ f{Xs)ds'^giXT] 



As in the proof of Theorem 15.11 we will use the parametrization in remaining time 
r = T — t to maturity. Let 

p(x,r)=e-^(-)^o^'"+''\ 

where the positive constant D satisfies @ and ([5]) for both models, and the constant 
Ko is large so that Kq > (e^^ -1)/D + Ke^^. For e > 0, define 

V%x,t) = V{x,t) + ee^^^^il + x^p{x,t)). 

It is straightforward to check that the constant M can be chosen so large that 

(33) M(l + x'^p) > -x^drP + adl{l + x^p) + (3dx{l + x^p) - /(I + x^p) 
at all points (x,r) G M x [0,T]. Now, let 

r := {(x,r) G M X [0,r] : V{x,t) < V^x,t)}, 

and suppose that T 7^ 0. Since Vp^^ is bounded by Lemma [3.21 and since V^p~^ grows 
like x^ for large |x|, there exists R > such that F C [—R, R) x [0, T]. Therefore, F is 
bounded and F is compact. Define 

To := inf{r : (x,r) G F for some x}. 

By compactness, there exists xq G M such that (xq, tq) G F. Since — V is continuous, 
we have y^(xo,To) = y(xo,To). Therefore, y^(x,0) > V{x,0) implies that tq > 0. By 
the definition of tq we have V'^{xq,t) — V{xo, r) > for < r < tq, so 

(34) dr{V^ -V)<0 

at (xo, To). On the other hand, V and V satisfy the parabolic equations Vr = C-^V and 
Vr = C-^V, respectively, where 

£/ = adl +(3dx-f and C = ad^ + pd^ - f. 

Consequently, at the point (xo,ro) we have 

dr{V'-V) = L^V - C^V + ee^^M{l + x'^p) 

+ee^'^ {x^drP -adl{l + x^p) - /39^(1 + x'^p) + /(I + x^p)) 

> aVx% + pV^ - - {aVxx + PVx - fV) 

by (|33p . The function x 1-^ V^{x,tq) — V{x,tq) attains its minimum at x = xq. Thus 
we have V^ = V, V^ = Vx and Vxx > Vxx at the point (xq, tq). Since at least one of the 
two models is convexity preserving, Vxx > at this point, so aVxx > aVxx- Therefore, 

at (xo,ro), 

dr{V' -V)> aV^x - aVxx + [[i-mx> 0, 
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where we also have used (3 < (3 and 14: < 0. This contradicts ()34p . Consequently, T is 
empty, so >¥ everywhere. Letting e \ it follows that V >¥ . 

Finally, an application of the monotone convergence theorem (as / approaches x) 
yields that U >U. □ 



7. Bond options 

In this section we apply the results of Sections [5] and [6] to study convexity and 
monotonicity properties of bond option prices. We consider a European call option 
with time of maturity Ti and strike price ivT > on a bond maturing at T2, where 
naturally T2 > 7i. We denote the price of this option at time t < Ti by C{x,t;Ti,T2) 
with X denoting the short rate. Thus 



C{x,t;Ti,T2) = E^, 



exp I - j^""' X, ds^ {uiXr, , Ti) - K)' 



where U{x, t) is the value function of a T2-bond. Assuming that the short rate dynamics 
is given by equation ([3]), we then have the following result. 

Theorem 7.1. Assume that I3xx{x,t) < 2 (in the sense of distributions). Then the 
bond call option price C{x, t; Ti,T2) defined above is convex in x at all times t £ [0, Ti]. 

Proof. We know from Theorem 15. II (applied with g = 1) that the price u of the T2-bond 
is convex in x. Note that an increasing function of the bond price is decreasing as a 
function of the short rate. Thus, since the pay-off function of the call is convex and 
increasing, it follows that C{x,Ti;Ti,T2) = {u{x,Ti) — K)~^ is convex. Moreover, since 
u is decreasing in x so is C(x, Ti; Ti, T2). Since there exists a constant M > such 
that 

{u{x,Ti)-K)+ < u{x,Ti) < Mmax{e^^,l}, 
see Corollary 13. 3| the result follows from another application of Theorem 15.11 □ 

We also have a related monotonicity result. To formulate this, let X and X be two 
diffusion processes satisfying 

dXt = p{Xt,t)dt + aiXt,t) dBt 

and 

dXt = P{Xt,t)dt + d{Xt,t)dBt. 

Let C and C be the corresponding call option prices written on the T2-bond prices u 
and u, respectively. Then we have the following result. 

Theorem 7.2. Assume that P{x,t) < P{x,t) and \a{x,t)\ > |fT(x,t)| for all x and t. 
Also assume that either Pxx < 2 (in the distributional sense) at all points or f3xx ^ 2 
at all points. Then 

Cix,t;n,T2)<C{x,t;TuT2) 

for t < Ti. 

Proof. First note that 

(35) g{x) := {u{x, Ti) - K)+ > {u{x, Ti) - K)+ =: g{x) 
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since u > u by Theorem 16. li If I3xx < 2 at all points, then g is decreasing and convex. 
It thus follows from Theorem 16.11 and (I35p that 



C{x,t;TuT2) 



> Ex,t 

> Ex,t 



exp 
exp 
exp 



Xsds}g{XT,) 



Xsds}g{XT,: 



Ti 



Xsds}~g{XT,: 



C{x,t;Ti,T2) 



A similar argument can be applied if instead Pxx < 2 at all points. This finishes the 
proof. □ 

Remark As noted in Sectional all models in Table [T] satisfy (3xx ^ 2. Thus all those 
models have bond call option prices which are convex, decreasing in the drift and 
increasing in the volatility. 



8. LOG-CONVEXITY OF BOND PRICES 

In this section we study convexity properties of the logarithms of bond prices. Re- 
call that a non-negative function u is said to be log-convex if u{\xi + (1 — X)x2) < 
u{Xxi)u{{l — X)x2) for all A G (0,1) and xi, X2- If u is strictly positive, then log- 
convexity is equivalent to the function x ^ Inn(x) being convex. 

For simplicity, we only deal with log-convexity of bond prices, so we assume that the 
pay-off function g = 1. Thus 



(36) 



W{x,t) 



,/(x){e^(^-')-l)/D 



E^ 



x,t 



exp 



T 



f{Xs) ds 



Recall from Lemma 13.21 that the function W is bounded on M x [0,T]. We first show 
that if the drift (3 is such that (jl5p is satisfied, then the function W is also bounded 
away from zero. 

Lemma 8.1. Assume that g = I, and that Hypothesis \ 3.1\ and (|15p hold. Then W 
defined in (j36p is also bounded away from 0, i.e. there exists a constant C > such 
that 



1 

C 



< W{x,t) < C 



on M X [0,T]. 



Proof. Using Jensen's inequality, applied to the exponential function and the expecta- 
tion, yields 

exp| re^(^-*)/(x)-/(X,)d. 



W{x,t) = Ex,t 

> exp <! Ex,t 



T 



M^-^)f{x)-f{Xs)ds 



exp|-^^e^("-*)£[n] dsj 



where again Y{s) = e f{Xs) — f{x) with initial condition Yj = 0. Note that the 

condition (jlSp yields that the drift /? of y is bounded also from above by some constant 
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C, SO E [Ys] < C{s - t) for s E [t, T]. Consequently, 

W{x, t) > exp 1^ -e^("-*)C7(s - t) dsj > expj-Ce^^T^} > 

for all t £ [0, T], which finishes the proof of the lemma. □ 

Theorem 8.2. Assume that a is spatially convex and (3 is spatially concave. Then the 
bond price u{x,t) is log-convex in the spot rate x. 

Remark It follows that all the models in Tabled] have log-convex bond prices. 

Proof. We need to show that lnu{x,t) is spatially convex. Similar to the proof of 
Theorem 15. 1^ we approximate the function In u by the function In 1/ = In V-^ for some 
function / satisfying Hypothesis 13.11 We first assume that Hypothesis 13.41 holds. In 
addition to the bounds on the derivatives in Hypothesis 13.41 we also assume that 

(37) \d^a{x,t)\ <C{l + x)^-'' A: = 0,1, 2 
for positive x. 

It is straightforward to check that the function P := InV-^ satisfies the non-linear 
equation Pr = CP, where 

CP = aP,, + + pp^ _ /, 

with initial condition P{x,0) = (here we have again used the parametrization in 
terms of time r to maturity rather than the physical time t). Moreover, there exists a 
constant C such that 

(38) \P.Ax,r)\<-^— 

for all (x,r) G M x [0,T]. Indeed, for large |x| we have fxx{x) = 0, so 

WWxx - wl 



Pxx = {lnV)xx = (Inl^ - fh)xx = (InW), 



which decays like according to Proposition 13.51 and Lemma l8.1i 

Let g : M — > M be a decreasing and convex function with strictly positive second 
derivative such that 

-x — (In \x\)'^ if X < —K 
-(Inx)^ if x > K 

for some positive constant K. We now claim that there exists a positive constant M 
such that 

(39) Mfe > d^iaq^x + PQx + ^aP^Qx) + \dl{aql)\ 

at all points {x,t) G R x [0,T]. To see this, note that qxx behaves like |rc|~^ln|x| for 
large \x\. The estimates 

Or . C , , , C 



\Px + /x(e^" - l)/D\ < ——- and 



1 + l^l 1 + I^P 

can be obtained in the same way as (j38p was derived. Using these estimates. Hypoth- 
esis [33] and (j37p . it is straightforward to check that all terms in 

dliaqxx + PQx + '^aPxQx) + \dl{aql)\ 

decay at least like |x|~^ln|x|. Consequently we can choose M large so that ([39]) holds. 
Next, for e £ (0,6"*^'^), define 

P^(x, r) := P{x, r) + ee^^q{x). 
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Since Pr = CP, straightforward calculations yield that 

mdl{P^ - CP') = ee^'^Mq^^ - ee^^dl{aq^^ + ee^'^aql + pq^ + 2aP,q,) 

> ee*^^ (Mfc - dliaq,, + Pq, + 2aP,q,) - \dliaql)\) > 0, 

where we used ee^^^ < 1 and ([39]). 

It follows from ([38]) and the fact that q^x behaves like x~'^ In for large |x| that the 
set 

r:={{x,r)GRx[0,T]:P^,{x,T) <0} 
is bounded. Thus, if F is non-empty, then there exists a point {xq, tq) G T such that 

tq = inf{r : (x, r) G F for some x £ M}. 
Since P^^{x,0) > we have tq > 0. Consequently, at the point (xo,ro) we have 

dlP^ = drPL < 0. 

Moreover, 

dUCP') = dl{aP^, + a{P^f + pP^-f) 

~ ^^xxxx ~l~ (2cix + 0)Pxxx ~l" {(^xx + '2f3x)Pxx ~^ Oixx{Px) 

-\-AaxPxPxx ~l~ '^'^■^x^xxx ~^ '^'^i-^xx) ~^ PxxPx ~ fxx 
> PxxP^>0, 

where the first inequality is due to Pxxxx ^ ^xxx = ^xx = at {xq,to), the assumption 
o^xx ^ and the concavity of /, and the second inequality follows from < and the 
assumption Pxx < 0. Thus 

(41) dl{P^-CP')<0 

at {xq,tq). But this is a contradiction to (j4U|) . which shows that the set F is empty. 
Thus is convex. Letting e tend to we find that also P is convex at all times 
rG[0,r]. 

Now, if Hypothesis 13.41 and (|37p are not satisfied, then we can approximate a and 
P with smooth coefficients as in the proof of Theorem 15. 1[ Using Proposition 14.11 and 
then letting / — > x it is straightforward to check that also Inii is convex, which finishes 
the proof. □ 

9. LOG-CONCAVITY OF BOND PRICES 

In this section we discuss concavity properties of the logarithm F = lnu{x,t) of the 
bond price. Note that F satisfies the non-linear parabolic equation 

(42) Fr = aFxx + aF^ + pFx - x 

with initial condition F{x, 0) = 0. In order to find the appropriate condition for 
preservation of log-concavity, we first present a heuristic argument similar to the one 
presented in the beginning of Section [5j 

Assume that {xq,tq) is a first point where concavity is almost lost, i.e. x F{x,t) 
is concave for all t ^ tq, and Fxx 

(xo,ro) = 0. Differentiating twice gives 
drFxx = dlFr = dl{aFxx + aF^ + pFx - x) 

— OiPxxxx ~l~ (2o^x + P)Fxxx ~l~ {oixx ~l~ ^Px^Fxx ~\~ PxxPx 
~\~OlxxPx ~^ ^OlxPxPxx ~\~ 2(y.{FxFxxx ~l~ -^xx)' 

Since Fxxxx < Fxxx = Fxx = at the point (xo,ro) we get 

drPxx ^ PxxPx ~l~ OlxxPx 



20 



ERIK EKSTROM AND JOHAN TYSK 



at that point. Thus, since < 0, we see that a sufficient condition for preservation 
of log-concavity appears to be that (3 is convex and a concave. Since a = a"^ /2 is 
non-negative, however, our convention that the model is specified on the whole real 
line is no longer convenient. Indeed, specifying a to be for negative short rates is 
not compatible with a concave infinitesimal variance a. The only possibility to have a 
concave on the whole real line is to require it to be constant in x. Such models can be 
shown to preserve log-concavity using the same methods as in Section [HI 

Theorem 9.1. Assume that a is only time- dependent, i.e. a{x,t) = ^{t) for some 
function 7. Also assume that is convex in x for each fixed time t. Then the T-bond 
price u{x,t) is log-concave in x at every fixed time t £ [0,T]. 

Proof. The proof follows along the same lines as Theorem 18.21 with some minor modi- 
fications. The function / of Hypothesis 13.11 needs to be replaced by a convex function 
which equals x for positive x and is constant for x < —K, where K is some positive 
constant. With this new /, Lemma l8. II and Proposition 13.51 remain valid if we modify 
P to be linear for x large positive. □ 

The conditions in Theorem l9.1l are only satisfied by the models V and HW in Tabled! 
To investigate the remaining models, we are forced to leave the tractable setting of 
diffusions defined on the whole real line and instead consider models specified on a 
half-line. Such models, however, typically lead to partial differential equations with 
degenerate coefficients. 

For simplicity, we assume that X is defined on the positive real axis [0, 00). We will 
also assume that a(0, t) = and /3(0, t) > 0. Note that under these conditions, no 
boundary behavior of the process at x = needs to be specified. Let 

(43) t/;(x,t) = e'^W^u(x,t), 

where h{t) = (e^('^-*)-l)/Z). By arguing as in the proofs of Lemma l3.2l and Lemma fH.H 
it can be shown that if jS—Dx and a are bounded for large x, then there exists a positive 
constant C such that 

(44) C-^<w{x,t)<C 

for all (x,t) G [0, 00) x [0, T]. To apply the techniques used in previous sections, we 
also need estimates of the derivatives of the function w. We have the following result. 

Lemma 9.2. Assume that X is specified on the positive real axis with q(0, f) = 0, 
/3(0, i) > 0, a{x,t) > for x G (0,cxd), that a and j3 are smooth and that a and 
(3 — Dx are constant in x for large x. Then there exists a constant C > such that for 
/c = 0, 1, 2, ... we have 

(45) |9^fc^(^,t)| <_il_ 
for {x,t) G (0,00) X [0,r]. 

Proof. As noted above, the case k = can be proven analogously to Lemma 13.21 It is 
straightforward to check that w := w — 1 is the bounded classical solution to 

(A(^\ ( wt + oiWxx + +7u;-F7 = t <T 

^ ' \ w = Q t = T, 

where 

$ = P - 2ah and 7 = ah"^ + {Dx - (5)h. 

Note that 

(47) \dli3{x,t)\<C{l + xf-^ and |5^7(x, i)| < C(l + x)-'= 
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for some constant C. By stochastic representation, 

l-T 

w{x, t) = E 



^ exp|^%(y,-'*,r)dr|7(y/'*)(i5 



where Y^'^ is a diffusion given by 

dyy^^ = /3(y/'*, s) + a(yf s) dB, 



and the indices indicate that X^' = x. As in the proof of Lemma 13.51 the bounds in 
(|47|) and Gronwall's lemma can be apphed to prove that Wx is bounded and satisfies 
Wx{x,T) = 0. Thus Wx is a bounded solution to the equation obtained by differenting 
the equation (|46|) . and using the maximum principle the estimate (|45p can be established 
for k = 1. The rest of the proof follows inductively by treating the differentiated 
equation as above. □ 

Remark It follows from Lemma 19.21 that the derivatives dxdj-u, k + 2l < A, are contin- 
uous up to the spatial boundary x = 0. 

Theorem 9.3. Assume that X is specified on the positive real axis with a{0, t) = 
and /3(0, t) > 0. Also assume that a is concave in x and (3 is convex in x at any fixed 
time t G [0,T]. Then the bond price u{x,t) is log-concave in x. 

Proof. We will assume that a and (3 satisfy the conditions in Lemma 19.21 We also 
assume that there exists a constant tj > such that a{x,t) = j{t)x and Pxx{x,t) = 
for X < r] and for some function 'y{t) > r]. The general case follows by approximation. 

Define the function F{x,t) = \nu{x,T), where n is as in ([T]). As above, F satisfies 
the non-linear equation 

Fr=£F 

where 

CF = aFxx + aF^ + /3Fx - x. 
It follows from ()44|) and Lemma 19.21 that 

\Fx{x,t) + h{t)\ <C{l + x)-\ 

\Fxx{x,t)\ < Cil + x)'^ 

and 

\Fxxxix,t)\ <C{l + x)-\ 
Let q : (0, oo) ^ M be a smooth, increasing and concave function with strictly negative 
second derivative such that 

X — x'^ ifx<l/C 
(Inx)^ if X > C 

for some constant C > 0. We claim that there is a constant M > so large that 
(48) Mqxx < dliaqxx + I3qx + ^aFxqx) - \dl{aql)\ 

at all points (x,t) G [0,oo) x [0,T]. Indeed, the right hand side is bounded for small 
X, and for large values of x the right hand side decays at least as fast as x~^\n.x. 
Consequently, M can be chosen so that (148 p holds at all points. 
Now, for £ E (0,e-*^^), define 

F^(x, r) = r) + ee^^q{x). 

Then 

mdl{F^-CF') = ee^'^Mqxx-dl{aqxx+ee^'^aql+Pqx + 2aFxqx)) 

< ee*^" {Mqxx - dl{aqxx + PQx + 2aF^g^) + \d^x(.aql)\) < 
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according to ([ISj) . Next, define the set 

r = {{x,t)eRx[0,T]:F^,{x,t)<0}. 

Since F^x decays at least like we have that F C [0, R) x [0, T] for some R, so F is 
compact. Let {xo,tq) be a point such that 

To = inf{r > : (x, r) G F for some x € [0, oo)}. 

If xq > 0, then arguing as before the inequality (j4ip . it is straightforward to check that 

dl{F^-CF^)>0 

at (xo,To), which contradicts (09]). Therefore, assume that xq = 0. Then, at the point 
{xq, To) we have 

(50) d^F^ = drF^, > 
and 

(51) dliCF') = dl{aF^, + a{F^f + (3F^-x) 

= O^xx^xx "I" '^'^x^xxx ~^ ^-^xxxx "I" ^xx^xx "I" '^'^xFx-^xx 

-\-2a{F^F^^^ + (F^^) ) + Pxx^x ~l~ '^PxF'xx ~^ l^-^xxx — 0' 

where we have used a = 0, /3 > 0, Fa;^; = and F^;!;!; < at the point (xo,to). The 
inequalities ()49p . (I5ip and ()50p form a contradiction. This shows that the set F is 
empty, so F^ is convex at all times. Letting e — > 0, it follows that also F is convex in 

X. □ 

Remark Prom the results in Sections [8] and [9] we find that if a and P are both concave 
and convex, i.e. linear, then both log-convexity and log-concavity are preserved. In 
our terminology such models thus give rise to log-affine bond prices. Of course, these 
models are usually referred to as having an affine term structure and they play an 
important role in interest rate theory. Our sufficient conditions for the existence of an 
affine term structure are well-known, see for instance Chapter 17 of [7]. However, the 
results of Sections [8] and [9] offer a background to these seemingly ad hoc conditions. 

10. Conclusions 

In this paper we have conducted a study of convexity of solutions to the term struc- 
ture equation. We show that if the drift P satisfies Pxx < 2, then the bond prices 
are convex in the current short rate, increasing in the volatility of the short rate and 
decreasing in the drift. Similar results hold for call options written on bond prices. For 
models with regular coefficients, the condition Pxx < 2 is also a necessary condition for 
preservation of convexity. We also have a general comparison theorem: if a model has 
smaller drift and larger volatility than another model, and at least one of them has a 
drift satisfying the condition above, then the first model has the larger bond prices. 
For bond call options the analogous result holds. 

We also study convexity properties of the logarithm of a bond price corresponding 
to the relative sensitivity of the bond price to changes in the short rate. This relative 
sensitivity is often described by the duration, i.e. the negative of the derivative of the 
logarithm. We show that if the drift /? is concave and the square cr^ of the volatility 
is convex, then bond prices are log-convex (a decreasing duration). Similarly, if /3 is 
convex and <t^ is concave, then bond prices are log-concave (an increasing duration). 
We also note that if we demand that the price is log-convex and log-concave, we recover 
the well-known sufficient conditions for a model to admit an affine term structure. 

Our findings for some commonly used models are summarized in Table [2] below. 



CONVEXITY THEORY FOR THE TERM STRUCTURE EQUATION 23 



Model 


Dynamics 


AB 


AO 


C 


LCV 


LCC 


V 


dX = k{9 - X) dt + a dB 


Yes 


Yes 


Yes 


Yes 


Yes 


CIR 


dX = k{e - X)dt + a^/XdB 


Yes 


Yes 


Yes 


Yes 


Yes 


D 


dX = hXdt + aX dB 


Yes 


No 


Yes 


Yes 


No 


EV 


dX = X{ri - alnX) dt + aX dB 


No 


No 


Yes 


Yes 


No 


HW 


dX = k{et- X)dt + adB 


Yes 


Yes 


Yes 


Yes 


Yes 


BK 


dX = X(r]t-a In X) dt + aX dB 


No 


No 


Yes 


Yes 


No 


MM 


dX-X{rit (A ^)lnX)dt + aXdB 


No 


No 


Yes 


Yes 


No 


Table 2. In the three last columns it is indicated whic 


1 moc 


els pre- 



serve convexity for option prices and bond call options, log-convexity of 
bond prices and log-concavity of bond prices, respectively. 
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